function output = ri(r,w,N,i,t)

hix_a = -w*w*r*cos(w*t + 2*pi*(i-1)/N);
hiy_a = -w*w*r*sin(w*t + 2*pi*(i-1)/N);

output = [hix_a;hiy_a];
